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(N : ABSTRACT 
> 

\Q I In this paper we compare the predictions of a detailed multi-zone chemical 

I evolution model for elliptical galaxies with the very recent observations of the 

O ' galaxy NGC 4697. In particular, the model allows for an initial gas infall and a 



subsequent galactic wind; it takes into account detailed nucleosynthesis prescrip- 



O ■ tions of both type II and la supernovae and reproduces the main photo-chemical 

' properties of normal ellipticals. As a consequence of the earlier development 

I ' of the wind in the outer regions with respect to the inner ones, we predict an 

?-j ■ increase of the mean stellar [< Mg/Fe >] ratio with radius, in very good agree- 

^ ■ ment with the data for NGC4697. This finding strongly supports the proposed 

outside-in formation scenario for ellipticals. 
^ In order to compare theory and observations it is necessary to convert either 

■ the predicted abundances into line strenght indices or the latter ones into abun- 

dances. We show that, in spite of the good agreement found for the [< Mg/Fe >] 
ratio, the predicted slope of the mass-weighted metallicity (Z) gradient does not 
reproduce the one derived from observations, once a calibration to convert in- 
dices into abundances is applied. This is explained as the consequence of the 
different behaviour with metallicity of the line-strength indices as predicted by a 
Single Stellar Population (SSP) and those derived by averaging over a Composite 
Stellar Population (CSP). In fact, in the current literature, the conversion indices- 
abundances is often obtained by assuming that ellipticals are SSPs, whereas they 
are CSPs at each radius. In order to better address this issue, we calculate 
the theoretical "G-dwarf " distributions of stars as functions of both metallicity 
([Z/H]) and [Fe/H], showing that they are broad and asymmetric that a SSP 
cannot correctly mimick the mixture of stellar populations at any given radius. 
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We also compute the stellar distribution as a function of the [Mg/Fe] ratio, which 
is narrower and with a negligible skewness than those as functions of [Z/H] and 
[Fe/H] and hence can be better represented by a SSP with an abundance ratio 
given by the average [< Mg/Fe >] ratio. Moreover, we compute the luminosity 
distributions of stars for a typical elliptical galaxy as functions of [Z/H], [Fe/H] 
and [Mg/Fe] ratios. We find that these distributions differ from the "G-dwarf" 
distributions especially at large radii, except for the one as a function of [Mg/Fe]. 
Therefore, we conclude that in ellipticals the [Mg/Fe] ratio is the most rehable 
quantity to be compared with observations and is the best estimator of the star 
formation timescale at each radius. 

Subject headings: galaxies: abundances — galaxies: elliptical and lenticular, cD 
— galaxies: formation — galaxies: individual (NGC4697) 

1. Introduction 

Metallicity gradients are characteristic of the stellar populations inside eUiptical galaxies. 
Evidences come from the increase of line-strength indices (e.g. Carollo et al., 1993; Davies 
et al., 1993; Trager et al., 2000) and the reddening of the colours (e.g. Peletier et al. 1990) 
towards the centre of the galaxies. The study of such gradients provide insights into the 
mechanism of galaxy formation, particularly on the duration of the chemical enrichment 
process at each radius. Metallicity indices, in fact, contain information on the chemical 
composition and the age of the single stellar populations (SSP) inhabiting a given galactic 
zone. In particular, by comparing indices related mainly to Mg (e.g. Mg^, Mg2) to others 
representative of the Fe abundance (e.g. < Fe >), it is possible to derive the [Mg/Fe] 
abundance ratio, which is a very strong constraint for the formation timescale of a galaxy. 
In fact, the common interpretation of the a-element (O, Mg, Ca, Si) overabundance relative 
to Fe, and its decrease with increasing metallicity in the solar neighbourhood is due to the 
different origin of these elements (time-delay model, Greggio & Renzini, 1983, Matteucci & 
Greggio, 1986), being the former promptly released by type 11 supernovac (SNll) and the 
latter mainly produced by type la supernovac (SNla) on longer timescales. The time-delay 
model applies also to other objects and the [a/Fe] versus [Fe/H] relation depends strongly on 
the star formation history. For a very short and intense star burst, the [ct/Fe] ratios decrease 
at higher metallicity than in the solar vicinity. The contrary occurs for a very slow star 
formation rate like in irregular galaxies (see Matteucci, 2001). Pipino & Matteucci (2004, 
PM04) showed that a galaxy formation process in which the most massive objects form faster 
and more efficiently than the less massive ones can explain the photo-chemical properties of 
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ellipticals, in particular the increase of [Mg/Fe] ratio in stars with galactic mass (e.g. Faber 
et al., 1992; CaroUo et al., 1993; Davies et al., 1993; Worthey et al., 1992). Moreover, from 
a preliminary analysis of metallicity and colour gradients, PM04 suggested that a single 
galaxy should form outside-in, namely the outermost regions form earlier and faster with 
respect to the central parts (see also Martinelli et al. 1998). A natural consequence of this 
model and of the time-delay between the production of Fe and that of Mg is that the mean 
[Mg/Fe] abundance ratio in the stars should increase with radius. We stress that the radial 
variation of the [Mg/Fe] abundance ratio has not been previously constrained in the study 
presented by PM04. On the other hand, an inside-out formation of ellipticals, as suggested 
to explain the abundance gradients in the Milky Way (Matteucci & Frangois 1989, Hon, 
Prantzos & Boissier, 2000, Chiappini et al. 2001), would produce a decrease of the [Mg/Fe] 
with radius, since the outermost regions would evolve slower than the inner ones. Therefore, 
it is very important to study the gradient of the [Mg/Fe] ratio, since this can impose strong 
constraints on the mechanism of galaxy formation. 

In this paper we compare PM04's best model results with the very recent observations 
for the galaxy NGC 4697 (Mendez et al. 2005, which also provide for the first time abun- 
dances measured from planetary nebulae (PNe) in an elliptical galaxy). The high quality of 
the optical spectrum gives us the opportunity to study in great detail the radial variation 
of the [Mg/Fe]. In order to compare the observations of NGC 4697 with the PM04 model 
predictions, we need to adopt a suitable calibration relation to convert our predicted abun- 
dances into line-strength indices. A preliminary study of the effect of the adopted calibration 
on the predicted radial gradients has been included in PM04. Here we present a deeper and 
self-consistent study on this topic, by using the same SSPs adopted by Mendez et al. (2005) 
to convert their observed indices into abundances. In particular, we show that the suggested 
outside-in formation process is supported by Mendez et al.'s (2005) data, independently from 
the uncertainties in the adopted calibrations. Moreover, we compare the radial gradients of 
abundances in stars as predicted on the basis of either mass-weighted or luminosity-weighted 
averages with observations. We show that the discrepancies arising from this comparison 
should be ascribed to the fact that a galaxy cannot be described as a unique SSP and that, 
at a given galactic radius, several SSPs contribute to the integrated spectrum. The paper 
is organized as follows: in Section 2 we briefly describe the adopted theoretical model; in 
Section 3 the predictions are compared with observations and in Section 4 some conclusions 
are drawn. 
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2. The model 

The chemical code adopted here is described in full detail in PM04, where we address 
the reader for more details. In particular, being NGC 4697 a IO^^Mq galaxy, we refer to 
PM04's Model lib for the same mass. This model is characterized by: Salpeter (1955) IMF, 
Thielemann et al. (1996) yields for massive stars, Nomoto et al. (1997) yields for type la 
SNe and van den Hoek & Groenewegen (1997) yields for low- and intermediate- mass stars 
(the case with tjagb varying with metallicity) . 

The model assumes that the galaxy assembles by merging of gaseous lumps (infall) on 
a short timescale and suffers a strong star burst which injects into the interstellar medium 
a large amount of energy able to trigger a galactic wind, occurring at different times at 
different radii. After the development of the wind, the star formation is assumed to stop 
and the galaxy evolves passively with continuous mass loss. We recall here that the assumed 
star formation efficiency is u = lOGyr"^, while the infall timescale is r = 0.4 Gyr in the 
galactic core and r = 0.01 Gyr at IRejf, respectively. These values were chosen by PM04 in 
order to reproduce the majority of the chemical and photometric properties of ellipticals. In 
particular, the radial decrease of the infall timescale is needed to reproduce at the same time 
both the average metaUicity and colour gradients of the observed sample of eUipticals. Only 
the predictions on the radial variation of the [Mg/Fe] ratio in stars have not been compared 
to the observations. We notice the all the PM04 subcases predict an a enhancement which 
increases with radius, although in different amounts. This is a consequence of the galactic 
wind developing earlier where the potential energy is lower, as already shown by Martinelli et 
al. (1998) in models without infall. Therefore, it depends only on the time delay and not upon 
specific model parameters. We recall that, at a given radius, both real and model galaxies 
are made of a Composite Stellar Population (CSP), namely a mixture of several SSPs, 
differing in age and chemical composition according to the galactic chemical enrichment 
history, weighted with the star formation rate. On the other hand, the line-strength indices 
are usually tabulated only for SSPs as functions of their age, metallicity and (possibly) a- 
enhancement. Therefore, in order to convert the predicted abundances for a CSP into indices 
(especially in the case of short burst of star formation), it is typically assumed that a SSP 
with a mean metallicity is representative of the whole galaxy. In other words, we made use 
of our knowledge about the chemical enrichment and the star formation histories to calculate 
either mass- or luminosity- weighted mean abundance ratios (i.e. [< a/Fc >], [< Fe/H > 
] and [< Z/H >]) for our CSPs, which are then used to derive the predicted indices for our 
model galaxies by selecting a SSP with the same values from the compilations available in 
the literature. In the following we will refer to mass averaged abundance ratios between the 
generic elements Xi and Xj with the symbol < Xi/Xj >. Following Pagel & Patchett (1975, 
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see also Matteucci, 2001): 

<Xi/Xj>=l rXi(S)/Xj(S)dS, (1) 
So Jo 

where Sq is the total mass of stars ever born contributing to the light at the present time. 
The V-luminosity weighted abundances will be denoted as < Xi/Xj >v. Following Arimoto 
& Yoshii (1987), we have: 

< Xi/Xj >v= J2 nk,i(Xi/Xj)iLv,k/ "k,iLv,k , (2) 

k,l k,l 

where Uk^i is the number of stars binned in the interval centered around (Xi/Xj)i with V-band 
luminosity Lv,k- Then we define the [< Xi/Xj >] ratios as log < Xi/Xj > — log(Xi/Xj)Q, 
taking the logarithm after the average evaluation (see Gibson, 1996). Generally the mass 
averaged [Fe/H] and [Z/H] are slightly larger than the luminosity averaged ones, except for 
large galaxies (see Yoshii & Arimoto, 1987, Matteucci et al., 1998). 



3. Results and discussion 
3.1. NGC 4697 

In Fig. 1 we show the predicted radial trends of both the [< Mg/Fe >] (solid line) 
and [< Mg/Fe >v] (dotted line) abundance ratios versus the observed one in NGC 4697. 
The latter is obtained by Mendez et al. (2005) by converting the line-strength indices into 
abundances by means of Thomas et al. (2003, TMB hereafter) SSPs. The agreement is 
remarkable, especially because we did not tune the input parameters (i.e. radius, mass) to 
exactly match NGC 4697. Both the observed increase of [Mg/Fe] with radius and the slight 
older age (~10 Gyr) of the outermost regions compared to the age of the inner core (~9 
Gyr) confirm PM04's model predictions, namely an outside-in formation process in which 
the central part of the galaxy form stars for a longer period compared to the most external 
regions. This can be explained in terms of galactic winds developing earlier where the local 
potential well is shallower (Martinelli et al. 1998). In Fig. 2 we compare our theoretical 
gradient for the total metallicity Z (solid line) with observations. The values for the inner 0.1 
Reff are in good agreement with observations, whereas our estimates for the outer regions are 
higher than the observed ones. This is because in the conversion from indices to metallicity 
it is adopted the approximation of one SSP and a CSP has a higher (mean) metallicity than 
a SSP with the same index (Greggio, 1997). This results from the metallicity distribution of 
stars and will be discussed further in the next section. A steeper slope (remarkably close to 
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Fig. 1. — PM04's model lib predictions for the mean mass-weighted [< a/Fe >] (sohd) and 
luminosity- weighted [< a/Fe >v] (dotted) abundance ratios in stars as a function of radius 
compared to the [a/Fe] derived for the galaxy NGC 4697 (Mendez et al. 2005, full squares) 
by means of TMB SSPs. 
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Fig. 2. — PM04's model lib predictions for the mean mass-weighted [< Z/H >] abundance 
ratio in stars (sohd hne) as a function of radius compared to the one derived for the galaxy 
NGC 4697 (Mendez et al. 2005, full squares) by means of TMB SSPs. For comparative 
purposes we also show the predictions for [< Fe/H >] (dotted line) and [< Z/H >v] (dashed 
line) . 
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the observed trend of decreasing [Z/H] with radius) is predicted in the case of [< Fe/H >] 
(Fig. 2, dotted hne). A better agreement with observation is also found for the predicted 
gradient of [< Z/H >v] (dashed hne, see Sec. 3.2). 

In order to better understand the origin of this discrepancy, here we focus on the proce- 
dure typically adopted to convert line-strength indices into real abundances. First of all, we 

run the procedure in the opposite sense, namely we use TMB's SSPs in order to transform 
our predicted [< Z/H >] and [< a/Fe >] into line-strength indices. The results for < Fe > 
are shown in Fig. 3 (we took also into account the radial age gradient suggested by Mendez 
et al. 2005). We obtain a good agreement for the < Fe > index, whereas the high predicted 
[< Z/H >] and [< Mg/Fe >] at large radii produce a positive (i.e. values increasing with ra- 
dius) gradient in the Mgt, index, at variance with observations. In fact, we obtain Mg^ — 4.1 
mag at 0.1 Reff, which is in good agreement with the observed value, but at 1 i?e// we have 
Mgi, — 4.6mag >> 3.7mag (observed). This discrepancy was already found by PM04 in the 
case of the Matteucci et al. (1998) a-enhanced cahbration. 

The reason for this disagreement is that the [< Z/H >] (which accounts for the ~ 60% 
of the predicted Mgi,, according to TMB) has an almost flat gradient (Fig. 2), whereas the 
[a/Fe] ratio (which contributes ~ 30% to the predicted Mgi,) increases with radius (Fig. 

1). The use of the luminosity- weigthed values helps in obtaining a negative gradient of Mgi, 
(namely Mgi, = 4.1 at 0.1 Reff and Mgi, = 4.0 at 1 Reff), but the slope is still too fiat 
with respect to the observations. As it will be shown in the next section, in fact, the use of 
the averages on the luminosity tends to increase the weight of the low metallicity SSPs with 
respect to the mass-weighted averages. On the other hand, we can still recover the gradient 
in the Fe indices (Fig. 3) because it is mainly driven by variations in the [< Fe/H >] alone, 
for which we predict a decrease with radius. 

3.2. CSP versus SSP 

Given the discrepancies encountered in the previous section, we further discuss what 
is behind the different behaviours of SSPs and CSPs. The main problem is that, in order 
to transform abundances into indices, one has to select a SSP which should be the closest 
approximation for the predicted CSP at a given radius. This is usually done by assuming 
that [< Z/H >] and [< a/Fe >] are representative of the galactic zone and, then, by 
interpolating among the values tabulated in a suitable grid. In the case of TMB's SSP, 
the relation [Fe/H] = [Z/H]-0.94[Q;/Fe] has been adopted, on the basis of a calibration on 
galactic globular clusters. This relation can be explained on theoretical grounds (Tantalo & 
Chiosi 2004) and holds also for the abundances in the gas of our model galaxy, but it is not 
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Fig. 3. — PM04's model lib predictions for the < Fe > line index (obtained by means of 
TMB's SSPs from mass-averaged abundances) as a function of radius compared to the one 
observed for the galaxy NGC 4697 (Mendez et al. 2005, full squares) . 
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necessarily true for a CSP. In particular, at 0.1 i?e// we predict [< Fe/H >] = 0.116 ^ [< 
Z/H >] - 0.94[< a /Ye >] = -0.18, whereas at 1 R^ff we predict [< Fe/H >] = -0.145 
[< Z/H >] — 0.94[< a/Fe >] = —0.52. Therefore, since the SSPs are calculated by means 
of fitting functions which depend on [Fe/H], this procedure introduces erroneously low Fe 
abundances. This is a general result and depends on the metallicity distribution of the 
stellar populations of an elliptical galaxy, as shown by Greggio (1997) for the indices in the 
nuclei of ellipticals. In particular, the calibration of line indices in terms of metallicity via 
the comparison with theoretical values for SSPs, is affected by a systematic effect if a CSP 
spanning a substantial metallicity range is present. 

In fact, from Fig. 4 it is clear that abundance ratios such as [a/Fe] have narrow and 
almost symmetric distributions, which encourages to take the average value as representative 
of the whole CSP. The robustness of the [a/Fe] ratios as constraints for the galactic forma- 
tion history is testified by the fact that [< a/Fe >] ~ [< a/Fe >v], having very similar 
distributions. On the other hand, the "G-dwarf" plots for [Z/H] and [Fe/H] exhibit a higher 
degree of asymmetry and have a broader distribution with respect to the [a/Fe] distribution. 
In particular, if we measure the asymmetry by means of the distribution skewness ^, we find 
that this parameter is much larger for the [Z/H] and [Fe/H] distributions than for the case 
of the [a/Fe] one by more than one order of magnitude. Moreover, the asymmetry increases 
going to large radii, especially for the [Z/H] distribution, in which the skewness increases 
by a factor of ~7 with respect to the inner regions. Therefore, it is not surprising that the 
[< Z/H >] value does not represent the galaxy at large radii, and hence, we stress that care 
should be taken when one wants to infer the real abundances of the stellar components for 
a galaxy by comparing the observed indices (related to a CSP) with the theoretical ones 
(predicted for a SSP). Only the comparison based on the [< a/Fe >] ratios seems to be 
robust. 

Another possible source of discrepancies is the fact that luminosity-weighted averages 
(which are more closely related to the observed indices) and mass-weighted averages (which 
represent the real distributions of the chemical elements in the stellar populations) might 
differ more in the most external zones of the galaxy. 

This was already shown in Fig. 2 (dashed line), where the radial gradient of [< Z/H >v] 
is in better agreement with the observed slope. In fact, in the luminosity-weighted averages, 
the stars with lower metallicity have a lower mass-to-light ratio in the V-band with respect to 
the high [Z/H] ones; therefore, the luminosity- weighted averages shift towards systematically 



^ Given a set of N measured values Xj with a mean /x, the skewness is Skew = jf J2jLii^^^^)^ ' where 
a is the distribution's standard deviation (Press et al., 1986). 
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lower values than the mass- weighted ones and the difference between the true [Fe/H] and 
the one obtained from the TMB calibration decreases. This is particularly true in the most 
external region (sec Fig. 4) , at variance with the typical assumption that [< Z/H >] ~ [< 
Z/H >v] for sufficiently old systems, which is still valid, however, for the galaxy as a whole 
(Matteucci et al., 1998, Gibson, 1996, see Sec. 2) and for the inner regions. 

Beyond one effective radius, PM04's best model predicts that both star formation and 
infall timescales are so short that a flattening of the radial gradients is expected. This result 
seems to be in agreement with the metallicity gradient inferred from PNe by Mendez et al. 
(2005), although they provide only lower limits. Wc hope that in the future more data from 
PNe (or other objects as globular clusters) will be available in order to better constrain the 
process of galaxy formation by means of new and independent tools. This could be a first 
step towards a better knowledge of the stellar abundances, as it happens for the stars in 
the Milky Way. In fact, the use of observations coming from single objects will avoid the 
use of the calibration relations with all the uncertainties linked to the population synthesis 
modelling. 

We conclude by stressing that NGC4697 is a well-observed single object, whose proper- 
ties have been compared to a model created to reproduce the average properties of ellipticals. 

Wc find that the agreement is remarkable. We notice that, recently, an increasing radial gra- 
dient for the [0/Fc] ratio has been observed by Proctor et al. (2005) in the galaxy NGC 821. 
The inferred radial trend is very similar to the one observed in NGC 4697 (see their Fig. 8). 
We avoided to show NGC 821 's observations, in order not to mix data obtained with different 
quality and calibrations. A larger compilation of [a/Fe] gradients is the one by Mehlert et al. 
(2003). They suggest that, at least for ellipticals in the Coma cluster, the average [Mg/Fe] 
gradient is zero (although objects with gradients as high as in NGC 4697 are also present). 
On the other hand, Ogando et al. (2005) presented observed massive ellipticals with strong 
gradients (up to A[Z/H]/ Alog{R/Reff) = —1) increasing with galactic mass. These slope 
are steeper than the typical values found by Mehlert et al. (2003) in their sample. For 
these galaxies we expect a very strong [a/Fe] gradient, but this should still be compared 
by observations. Massive galaxies which exhibit shallower gradients, instead, could have a 
prolonged formation history, perhaps affected by mergers (e.g. Kobayashi, 2004). 

Given the importance of the subject in the study of galaxy formation, and the issues 
arising when SSPs are taken as representative of CSPs, we hope that wider, more accurate 
and homogeneous samples of [Mg/Fe] gradients will be the goal of future observational 
campaigns. 

Finally, we notice that radial variations in the a enhancement alone can be reproduced 
in several ways. For example, a constant [< Mg/Fe >] along the galactic radius can be a 
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Fig. 4. — "G-dwarf" distributions for [Mg/Fe] (top row), [Z/H] (middle row) and [Fe/H] 
(bottom row) for the model galaxy in luminosity (solid line) and mass (dotted line). Left 
column: values at O.l-Rg//- Right column: values at 1-Re//- The plots are presented in the 
same scale in order to better appreciate the differences among the different distributions. 
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result of either a formation process which has a timescale independent from the galactocentric 
distance or by a radial change of the IMF (see e.g. Martinelh et al. 1998). In any case, it 
is the whole picture coming from the largest possible set of spectro-photometric observables 
which helps is determining the most likely formation scenario. 

4. Conclusions 

By means of the comparison between PM04's best model predictions and the recent and 
high quality data on NGC 4697 (Mcndcz ct al., 2005), wc arc able to derive some conclusions 
on the radial trend of abundances and abundance ratios in the stellar component of elliptical 
galaxies. Our main conclusions are: 

• PM04's best model prediction of increasing [< a/Fe >] ratio with radius is in very good 
agreement with the observed gradient in [a/Fe] of NGC 4697. This strongly suggests an 
outside-in galaxy formation scenario for elliptical galaxies that show strong gradients. 

• By comparing the radial trend of [< Z/H >] with the observed one, we notice a 
discrepancy which is due to the fact that a CSP behaves in a different way with respect 
to a SSP. In particular the predicted gradient of [< Z/H >] is flatter than the observed 
one at large radii. Therefore, this should be taken into account when estimates for the 
metallicity of a galaxy are derived from the simple comparison between the observed 
line-strength index and the prediction for a SSP, a method currently adopted in the 
literature. 

• Luminosity-weighted abundances might differ from mass-averaged ones at large radii 
even for old (i.e. ~ 10 Gyr old) stellar populations, at variance with what is commonly 
assumed. 

• Abundance ratios such as [Mg/Fc] arc less affected by the discrepancy between the 
SSPs and a CSP, since their distribution functions are narrower and more symmetric. 
Therefore, we stress the importance of such ratio as the most robust tool to estimate 
the duration of the galaxy formation process. 
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